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Abstract 

We develop an underlying relationship between the theory of rational approximations 
and that of isomonodromic deformations. We show that a certain duality in Hermite’s two 
approximation problems for functions leads to the Schlesinger transformations, i.e. trans¬ 
formations of a linear differential equation shifting its characteristic exponents by integers 
while keeping its monodromy invariant. Since approximants and remainders are described by 
block-Toeplitz determinants, one can clearly understand the determinantal structure in isomon¬ 
odromic deformations. We demonstrate our method in a certain family of Hamiltonian systems 
of isomonodromy type including the sixth Painleve equation and Gamier systems; particularly, 
we present their solutions written in terms of iterated hypergeometric integrals. An algorithm 
for constructing the Schlesinger transformations is also discussed through vector continued 
fractions. 
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Introduction 

Let L be an integer greater than one. For a given L-tuple of analytic functions (or formal power 
series) /o(w), /i(w),..., /l-i(w) and of nonnegative integers m = (mo, mi,..., m^-i), Hermite con¬ 
sidered the following two rational approximation problems. The first is to find L polynomials 

qr(w) (0</<L-l) 


of degree m, - 1 such that 


L-l 


1=0 


where \m\ = mg i.e., the left-hand side has a zero of order at least |»i| -1 at w = 0. The second 
is to find L polynomials 

vTiw) (0<i<L-l) 


of degree \m\ - m,- such that 

Each system of polynomials {q™} and {p("} generically turns out to be unique up to simultaneous 
multiplication by constants, as an elementary consequence of linear algebra. The above approxi¬ 
mation problems come from Hermite’s study on arithmetic properties of the exponential function 
and are called collectively the Hermite-Pade approximations', often, the former is referred to as the 
‘type r and the latter as the ‘type IT or as the simultaneous Fade approximation. Note that if L = 2 
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then both of them reduee to the (usual) Fade approximations. Although these two types of approx¬ 
imations were seemingly unrelated, Mahler diseovered that they were fundamentally eonneeted to 

eaeh other. Put e, = (0,..., 0, 1,0,..., 0) 6 for 0 < z < L - 1. 

Theorem 0.1 (Mahler’s duality). It holds that 
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with D being a diagonal constant matrix. Moreover, if every diagonal part and p'" is chosen 
to be a monic polynomial then D becomes the identity matrix. 

The aim of this paper is to develop an underlying relationship of Hermite’s two approxima¬ 
tions with the theory of linear differential equations in the complex domain, especially with that of 
isomonodromic deformations. Interestingly enough, Mahler’s duality plays a crucial role in con¬ 
structing a certain class of Schlesinger transformations, i.e. transformations of a linear differential 
equation shifting its characteristic exponents by integers while keeping its monodromy invariant. 

In Sect. [H we begin by introducing the two types of rational approximations for an L-tuple 
of functions, which are slightly modified (in order to fit the construction of Schlesinger transfor¬ 
mations) from the original Hermite-Pade and simultaneous Fade approximations. We then prove 
a variation of Mahler’s duality between them (see Theorem 11.31) . Applying the approximations 
for the solution of an L x L Fuchsian system of linear differential equations yields its Schlesinger 
transformation (see Theorem 11.61) : in fact, Mahler’s duality guarantees the absence of apparent 
singularities in the new Fuchsian system after the Schlesinger transformation. In Sect. we de¬ 
duce determinantal representations for the approximants and remainders from the approximation 
conditions (see Propositions 12. II and 12.31 and also RemarkIt should be noted that any key in¬ 
gredient here is written in terms of block-Toeplitz determinants. In Sect.|3l we present an algorithm 
for constructing the Schlesinger transformation via vector continued fraction expansions, which is 
a variation of the Jacobi-Perron algorithm (i.e. a higher dimensional analogue of the Euclidean 
algorithm). 

The last two sections are devoted to the study of isomonodromic deformations. In Sect. IH 
we first review the Schlesinger system of nonlinear differential equations, which governs isomon¬ 
odromic deformations of a Fuchsian system. Since a Schlesinger transformation preserves the 
monodromy of the Fuchsian system, it gives rise to a discrete symmetry of the corresponding 
Schlesinger system. We clarify, based on the above relationship with rational approximations, the 
determinantal structure in the general solutions of the Schlesinger systems. Next we concern a par¬ 
ticular family of the Schlesinger systems, which possesses a unified description as a polynomial 
Hamiltonian system denoted by PIl.n {L > 2, N > 1); it includes various noteworthy examples of 
isomonodromic deformations such as the sixth Painleve equation ('^ 2 , 1 ) and the Gamier system 
in N variables ("7^2 ,a?). In Sect. [H we demonstrate Schlesinger transformations on the previously 
known hypergeometric solution of ^l,n (see Il25l ): as a result, we obtain solutions of ^l,n writ¬ 
ten in terms of iterated hypergeometric integrals through Fubini’s theorem and the Vandermonde 
determinant (see Theorem 15. 3 1 and its sequel). 
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1 From Hermite’s two approximation problems to Schlesinger 
transformations 

In this section we show how rational approximations are useful for constructing Schlesinger trans¬ 
formations of linear differential equations. Fix an integer L > 2. We shall first introduce two 
different types of rational approximation problems for an L-tuple 

/o(w), /l(w), . . . , /l_i(w) 6 CM 

of formal power series, where we assume /o(0) 0 without loss of generality. 


1.1 Hermite-Pade approximation (of type I) 


Let n be a positive integer. Consider for each / (0 < i < L - 1) an L-tuple of polynomials 
Qf = (0 ^ j < L - 1 ) of degree at most n - I + 6ij, where 




ij - 


1 (i = j) 
0 (i j) 


is the Kronecker delta. Suppose the approximation condition 

2o^/o + • • ■ + Qf fi + wQflifi+\ +-1- wQf_\fL-\ = 0(w"^) 


(1.1) 


is fulfilled for each i. This condition amounts to a system of nL homogeneous linear equations for 
the nL+\ unknown coefficients of the polynomials Qf (0 < j < L-l). Under a generic condition 
for the power series /o, ... these polynomials are uniquely determined up to simultaneous 
multiplication by constants; see Sect. 12.11 We will be concerned with the row vector 




0<j<L-\ 


( 2 : 


Q , ■ ■ ■ ,Qj , 



Remark 1.1. By construction, the polynomial Qf^iw) has no constant term. Moreover, the degree 
of the diagonal part Qf{w) (0 < z < L - 1) turns out to be n exactly; see (11.41) in Sect. 11.31 


1.2 Simultaneous Fade approximation 


We treat another type of approximation problem for the same power series /o,..., /l-i. Consider 
for each 7 (0 < j < L - 1) an L-tuple of polynomials Pf = Pf\w) (0 < z < L - 1) of degree at 
most zz(L - 1) - 1 -I- 6ij which satisfies the following approximation conditions: 

• if 7 = 0 foPf^ - fiPf^ = Oiw"^) fori <z< L-l; 

• if 1 < 7 < L - 1 foPf - fiPf^ = for 1 < z < 7 - 1, (1.2) 

hPf - fiwPf^ = 0(w”^) for 7 < z < L - 1 . 


These conditions are interpreted as a system of nL(L - 1) homogeneous linear equations for the 
nL(L - 1) -I- 1 unknown coefficients of the polynomials Pf (0 < z < L - 1). Hence the column 
vector 



T 



0<(<L-1 



.,wPf„Pf, 



is generically unique up to multiplication by constants; see Sect. 12.21 
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Remark 1.2. Let a b.\{. is immediate from l//o 6 C[[w]] to verify that 

• if a,b<j 

• otherwise faP^l^ - ftP^a = 0 (w”^). 

1.3 Mahler’s duality in the two approximation problems 

There is an interesting connection between the two approximation problems (11.11) and (11.21) al¬ 
though they are seemingly unrelated. The following theorem is thought of as a variation of 
Mahler’s duality, see Theorem 10.II or ifTSll . We will give a proof of it because our setup is slightly 
different from the original; cf. ^ Theorem 8.1.2] and [|5l . 

Theorem 1.3. It holds that 


( 0 ) 


Q 


~(L-1) 

Q 




= w"^-D 


with D being a diagonal constant matrix. 

Proof. Let us first estimate the degree of the (z, y)-entry 


L-\ 

' ’ -i(0 n(y') 


^7 = 2 . Q'kP'k 

k=0 


of the left-hand side. The degree of each polynomial 2^'^ or P^J'' reads 


deg < H - 1 deg < n 

pU) p(J) pU) 

wPq ,... ,wF ■ , F . ,F^_^ . 


yo r>(i) 


,(0 


i (0 


deg < n{L — 1) 


deg < n(L - 1) - 1 


Hence we have 


deg^ Mij < 


nL (z < j) 
nL - 1 (z > j). 


(1.3) 


Next we shall estimate the multiplicity of M,y at w = 0 by means of the approximation condi¬ 
tions. Consider 


Z.-1 




(0 pij) 


k=Q 

= AQ'fPJ' + j; St APT- 

ki=0 
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(i) Case i < j (strictly upper triangular part) 


E 2“’/of = E +E i'’ 

ki:0 0<k< j j<k 

= Yu Qk + Xi Qf (fkPo'’ + Oiw'^'^)), using (D 

0<k<j j<k 

= ^ Qffk + since 2 ^'^ is divisible by w if / < (7 <)k. 

Therefore, 

foMij = Yj Qkfk + 

k 

= • 6 >(w”^) + 6 >(w”^+^), using (O) 

= 0{w''^^^), since is divisible by w if 0(< i) < j. 

(ii) Case i > j (lower triangular part) 

E SVo?i'’ = E +E Sffi'f’f 

ki^O 0<k< j j<k 

= Y Qk + 0 (w”^"i)) + ^ Qf (fkPf + 0 (w”'-)), using (O) 

0<k<j j<k 

= Po^YQkfk + 0(w'^^). 

ki=0 

Therefore, 

foMij = P^^^YQkfk + 0(.w^^) 

k 

= P^^ • 6 >(w"^) + 6 >(w”^), using ^ 

= 6 >(w"^). 

Noticing 1 //o 6 C[[w]] we verify 

( o(w”^+i) a < j) 

\ o(w”^) a > j). 

Combining this with (I1.3L we can conclude that (M,y) = ■ D. □ 

Consequently, the diagonal entry Mu coincides with the term of highest degree in (^^{w)Pf{w) 
and thus 

deg„ Qf = n , deg^ = n{L-\). (1.4) 

We henceforth normalize ^ and so that their diagonal parts and P^j^ become monic 
polynomials, i.e. 

and thereby D = I (the identity matrix). 
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Corollary 1.4. The polynomial matrix 


R{z) = z" 


r(0) 


Q iz-^) 


Q {z-^) 


■nLxL 


[z] 


(1.5) 


satisfies 


(i) R(zr^ = z 

(ii) det7?(z) = 1 


.n(L-l) 


~( 0 ) , , 
P\z-^),...,P V) 

/.e. 7?(z) 6 SL(L,C[z]). 


Proof. Theorem 11.31 shows (i) immediately. Then, it holds that R 6 GL(L, C[z]) and thus deti? 6 
C \ {0}. By definition, R takes the form 


1 • • • * 



+ 0{zf 


1 


( 1 . 6 ) 


namely, its eonstant term is an upper triangular matrix whose diagonal entries are all one. There¬ 
fore, we have detR =1. □ 


1.4 Schlesinger transformations 


Consider an Lx L Fuehsian system 


dT 

dz 



(A, : eonstant matrix) 


of linear ordinary differential equations with -l- 3 regular singularities 

S = {Mo = 1, Ml, . . . , Un, Un+I = 0, Un+2 = oo} c = C U {oo} 


(1.7) 


on the Riemann sphere. Let and An +2 = - ZZo' upper and lower triangular matriees, 
respeetively. Assume for simplieity there is no integer differenee among the characteristic expo¬ 
nents {eo,y}o<i<L-i at z = 0 (resp. {eoo,;}o<;<L-i at z = oo), i.e. the eigenvalues of the residue matrix 
A^+i (resp. An+ 2 )- Then we have a solution Y = y(z) of (11.71) normalized as 


VL 


1 ... * 


1 


+ 0(z) 


diag (z''"'^)o<;<L-i 


= 0(w). diag (w''”'Oo<;<L-i • C 


( 1 . 8 ) 


with w = l/z and C being an invertible eonstant matrix (the connection matrix between z = 0 and 
z = oo). Here 0(w) is a matrix funetion holomorphie at w = 0 (z = oo) and 0(0) is invertible and 
lower triangular, i.e. 

* 


0(w) = 


-I- 0{w). 


1 


















Remark 1.5. Many literatures adopt a different normalization sueh that the residue matrix at z = oo 
beeomes diagonal. Our present normalization treats the two points z = 0 and z = oo equally and it 
emerges naturally from the similarity reduetion of the UC hierarehy, whieh is a eontext of infinite¬ 
dimensional integrable systems; see [|^|^[26l. Furthermore, as elarified by Haraoka |l6l, this 
normalization is effective to find a ‘good’ coordinate of the space of Fuchsian systems having a 
given Riemann scheme. 

An analytic continuation along a loop on P' \ 5 based at some point zo induces a linear trans¬ 
formation of Y according to its multi-valuedness at the branch points S. We thus obtain an L- 
dimensional representation of the fundamental group 7ri(P^ \ S ; zo), which is called the monodromy 
of Y. A left multiplication Y Y = RY of a rational function matrix R = R{z) is said to be a 
Schlesinger transformation if the new equation 



satisfied by Y becomes the same form as the original (11.71) . Because R{z) is rational, Y and Y 
have the same monodromy though they have different characteristic exponents by integers. It 
is known that if we specify an admissible discrete change of the characteristic exponents then 
the corresponding rational function matrix R of the Schlesinger transformation is algebraically 
computable from 7; see |l9l[23. In fact, the construction problem of Schlesinger transformations 
is naturally related to rational approximation problems; see also Remark [L8l 

In this paper we focus on a class of Schlesinger transformations, which is of particular inter¬ 
est from the viewpoint of Hermite’s two approximation problems and also of vector continued 
fractions (see Sect. [3]). Note that, for a general Schlesinger transformation other than the present 
direction, though it can also be controlled by some rational approximation problems but it becomes 
much more complicated due to the absence of a duality like Mahler’s; e.g. R~^ seems not to have 
a concise determinantal representation. 

Let us define the L-tuple / = ^(/o,..., fi-i) of power series in w as the first column of 0(w), 
where 0(w) is the power series part of the solution Y of (11.71) near z = oo. Notice that /o(0) ^ 0 
certainly holds. Therefore, all the general arguments in Sects. [TTTI - ll.Sl are still valid for this specific 
case, and we are led to the Schlesinger transformation through the two approximation problems 
for /. Now we state the result. 

Theorem 1.6. The polynomial matrix R{z) given by (11.51) realizes the Schlesinger transformation 
Y Y = RY shifting the characteristic exponents at z = oo by (n{L - 1), -n ,..., -n). 

Proof It follows from the Hermite-Pade approximation condition (11.11) that Rf = By 

definition, R takes the form 




0 



R = 


* ^ 


-1- 0{w) 


1 


* * 

/ 


(1.9) 


Hence, if we write as 

R<I>(w) = 6(w) • diag ..., w“"). 
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then 6 becomes the same form as O. Also, we verify from (11.61) that R does not change the form 
of the power series expansion of Y near z = 0. 

On the other hand, the coefficient 


N+l 


A(z) 


V 


of the Fuchsian system (11.71) is transformed as 

-V , dR , 

A A = RAR-^ +—R-\ (1.10) 

dz 

If we remember both R and R~^ being polynomials (see Corollary 11.41) . then A turns out to be a 
rational function matrix having only simple poles at S as well as the original A. In this sense 
Mahler’s duality guarantees the absence of apparent singularities in the new equation dF/dz = AY 
satisfied by 7 = RY. □ 

Remark 1.7. In the rank two case (L = 2) a similar construction of Schlesinger transformations as 
Theorem 1 1.61 has been established in ifT^ based on (usual) Fade approximations. 

Remark 1. 8 . A series of pioneering works was done by D. Chudnovsky and G. Chudnovsky on the 
close connection between rational approximation problems and Riemann’s monodromy problem, 
involving (semi-classical) orthogonal polynomials; see [|3l|4l and references therein. The ‘Fade 
method’ recently proposed by Yamada [|2^ is a recipe for Lax formalism of Fainleve equations 
and, at the same time, for their special solutions, which is based on Fade approximations (or inter¬ 
polations) of elementary functions; interestingly enough, it is applicable also for various discrete 
analogues of Fainleve equations beyond the originals; see [rTl fTTlflgll . 

The essential idea of the above works could be exemplified by the following: let us consider a 
function ^(z) = z''(z-1)^(z-m)'^ witha-t-Z^-l-c = 0. The remainder p := P(p-Q = 0(z“”"^) of its Fade 
approximation then satisfies a second-order linear differential equation denoted by E, which may 
have an apparent singularity besides the four regular singularities S = {0,1, w, oo} c P'. However, 
the two functions ^ and p share the same multi-valuedness since they are rationally related; thus, 
the monodromy of E is obviously constant with respect to u. This fact leads to special solutions of 
the sixth Fainleve equation Fyn i-e. the isomonodromic deformation (cf. Sect.lH) of a second-order 
linear differential equation with four regular singularities. 

It is interesting to note that such an idea had been recognized implicitly by Laguerre (before 
the discovery of Fainleve equations); see [[T3ll and also ifTdl . 

Remark 1.9. The approximation conditions (11.11) and (11.21) can be interpreted as certain multi¬ 
orthogonality relations among the L-tuples of polynomials Q and respectively; i.e., these 
polynomials can constitute multi-orthogonal polynomial systems. In this paper, although we do 
not enter into details on such aspects, we present below the determinantal representations for them, 
which will crucially work in the last two sections. 
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2 Determinantal representations for approximation polynomi¬ 
als and remainders 

In this section we derive determinantal representations for the approximation polynomials Q^l\w) 
and P^j\w). We write the power series as 


fiiw) = ^ a‘jW^ 6 CM (2.1) 

j=o 

henceforth; note that the superscript i of a'j is just an index, not an exponent. Introduce the kx I 
rectangular Toeplitz matrix 

^jiktl) — \tt^j+m-n\l<m<k 
l<n<l 



1 

2 

1 

1 

a‘ 

a)_, ■ 


2 


q}- 

Uj 


k 

- 

^j+k-2 



for the sequence where a‘j = 0 if y < 0. It holds that 

(A^(m,n)), , = (Af„(n,m)) (2.3) 

by definition. 


2.1 Hermite-Pade polynomials 

We can calculate separately for each / (0 < / < L - 1). Therefore, for brevity, we shall express the 
coefficients of the approximation polynomial QJ\w) as 

n-l +djj 

e®(w)= ^ 

k=0 


with omitting the index i. The condition (11.11) implies vanishing of the coefficients of 1, w, ..., ^ 

in the left-hand side. Consequently, we have a system 



bo 


0 

' 


by 

= 

0 

>- 


br-y 


0 



of homogeneous linear equations for the nL + 1 unknowns 


(2.4) 


■ ■ ■ ,bj^n-\+6i,j) (0 < J < L - 1) 
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where 




(i) 


AQ(nL, n) 


A^\nL,n) A'^{nL,n + 1) A'll{nL,n) 


A^\\nL,n) 


i + 1 blocks 


L - / - 1 blocks 


The solution of (12.41) is unique up to multiplication by constants if and only if the rank of the 
nL X {nL +1) matrix equals nL (which we will always assume). 

Interestingly enough, we have the following determinantal representation of 

Proposition 2.1. It holds that 


QfM 


NQ' 


,('■) 


det 


jth block 


0 0 




(2.5) 


where are some normalizing constants. 

Proof. Consider 

Pi{w) = QqVo -I-+ Qf fi + -I-1- 

which is the remainder of the approximation condition (11.11) . Substituting (12.51) shows that 


Piiw) 


NQ' 


(0 


det 


0th block 


(/ — l)th block 


ith block 


/o, fow,..., /ow' 
A^(nL, n) 


n-l 


f.ufi-lW,...Ji-,W^-^ fJiW,...,fw” 


AL ^ (nL, n) 


(i + l)th block 
,2 


f+iw,f+iw ,...,f+:w" 
Af^ (nL, n) 


A’^inL, n + 1) 

{L— l)th block 


A^^^(nL, n) 


Therefore, if we put p,(w) = then the coefficients read 



1 


Al(nL, 


A}f\\,n) A[(l,n+\) At\(\,n) 

A‘f^(nL,n) AQ(nL,n + \) Af^(nL,n) 


A[:!(l,n) 

A^^^(nL,n) 


It is immediate from a property of determinants to verify p[ = 0 for any k less than nL; thus, we 
have Pi(w) = 0{w"^) indeed. □ 


We will normalize the polynomials so that its diagonal part 2f\w) becomes monic as well as 
in Sect. 11.31 Accordingly, the normalizing constant NQ^'^ should be 


NQ^'^ = det 


ith block 

0 oTT^Ti 0 


= (_1)«('+!) (Jet |^Ap(nL,n) ••• AQ(nL,n) A‘f^(nL,n) A^^^(nL,n)^. 
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Thus, the leading coefficient of the remainder is given by 

det [ A^Q{nL +l,n) A'~^{nLl,n) A'^^{nL + \,n + \) 

A‘*l{nL+l,n) A[_\\nL + l,n) 

the constant term of the polynomial Qf{w) (i 0) is given by 


e''’(0) = ^det 


1 


NQ(') 
(- 1 ) 


ith block 

0 iX^TTo 0 




— det|^AQ(nL,n) ••• A'^^{nL,n) A'_^{nL,n) ••• n)j 


NQ 

Remark 2.2. We here restrict ourselves to the case where /o(w) = 1. Let us introduce the block- 
Toeplitz determinant 


= det 


A^(m,n) ■■■ A\ (m,n) A'„_^{m,n) A„_i(m,R) 


i - 1 blocks L - i blocks 

of size m = n{L - 1) for each i (1 < z < L); e.g. We see in particular that 

AW 


pIl = (- 1 )" 




Ad)’ 

A(') 


These simple formulae will be used later in Sect. [51 


2.2 Simultaneous Fade polynomials 

Suppose /o(w) = 1 for simplicity. Or, equivalently, we may understand that we have renamed fi/fo 
(i 4^ 0) as fi. For a given power series F{w) = we employ the notation 

b 

VF{w)t = YjPk^'' 

k=a 

denoting its section between and if a < b. From now on, we set 

m = n{L - 1) 


as well as in Remark |2!^ 

First we shall construct the formulae for (0 < j < L - 1). 
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(i) Case j = 0 The approximation eondition (11.21) requires that 


[/cir'=o (i<,-<L-i) 

sinee Pf^ (i ^ 0) is a polynomial of degree at most m - 1. If we write 

m 

k=0 


then we find a system 


Al(n,m+ 1) 

Al(n,m+ 1 ) 

+ 1 ) 



( 2 . 6 ) 


of homogeneous linear equations for the m +1 unknowns bo,, b,n. The solution of (12.61) is unique 
up to multiplieation by eonstants if and only if the rank of the m x (m + 1 ) matrix in the left-hand 
side equals m. 


(ii) Case 1 < j < L - I Similarly, it follows from (11.21) that 

nm+rt -1 

0 (l<f<y-l) 
= 0 (/ = j) 


[f:P: 
[/jfo’l 


m 

m+n-2 

m 

im+n-2 


These amount to the simultaneous linear equation 


bo 


0 

' 

bi 





0 


bm-\ 





> m — 1 


A^~^{n,m) 
Ai{n - l,m) 


for the m unknown eoeffieients bo,..., b,„-\ of the polynomial 

m -1 

P^i'{w) = Y,bkwK 

k=0 

Proposition 2.3. The polynomials admit the following determinantal representations'. 


Pfiw) 


Np( 0 ) 


det 


l,w, ... ,w” 
-I- 1) 


\n,m+ 1 ) 
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and 


Pfiw) 


NP' 


ii) 


det 


l,w, 

A^(n, m) 

Ai^\n,m) 
Ai(n - l,m) 

for 1 < i < L - 1, where are some normalizing constants. 
Next the other (i 0) ean be written as follows: 


•if7 = 0 
•ifl< j < L - I 


piO) ^ 

pU) _ 


f. 

JOo 

or 

0 




m-l 


0 , 
m-l 


= w 

J 

P^j^ = w 


f.p'f 


m-l 


0 , 
m-2 


for 1 < / < j - I, 
for i = j, 

for 7 + 1 < z < L - 1. 

3(0/ 


We will ehoose the normalization so that the diagonal part P. \w) beeomes monie. Aeeordingly, 
we obtain 


NP(°^ = det 

0 ,...,0,1 
Ain(n,m+ 1) 

= (-l)'”det 

Ai„(n,m) 


A^-fn,m+ 1) 


Ainfn^m)^ 




(2.7) 


and 


NP^.'* = det 


a^ ,,a^ ai 

m-l’ m-2’ ’ 0 

A^(z7 , m) 


Ain \n,m) 
Am(zz - l,m) 

A,C-i («’"») 

A^:\(n,m) 


(_l)«(;-i)det 


Alfn, m) 

Aif\n,m) 

A;^_i(n,zw) 

Af„:'i(zr,zn)J 


( 2 . 8 ) 


for 1 < j < L - Here we have employed the notation of the bloek-Toeplitz determinant (see 
Remark [2^ in view of (12.31) . 

Remark 2.4. If z > y > 1, P\^\w) is divisible by w. By eonstruetion (see Corollary 1 1.4L we observe 
that takes the form (ef. (11.91) ) 


R-' = w- 


VL 


0 


+ 0(w) 
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3 From vector continued fractions to Schlesinger transforma¬ 
tions 

In this section we present an alternative construction of the same Schlesinger transformation (con¬ 
sidered in Sect. [B through an algorithm for expanding a vector-valued function into a vector con¬ 
tinued fraction. 


3.1 Algorithm for vector continued fraction expansion 

Let / = ^(/o,... be the L-tuple of formal power series (12.11) . For simplicity, we assume 
tentatively that fi(0) 0 for all i. We abbreviate the constant term aj, of fi(w) = 2/lo 

First we apply a left multiplication of a permutation matrix to /: 


0 1 


fo 


/l-1 

1 0 


fl 

— 

fo 

1 0 


fL-l_ 


1 

to 


Next we eliminate the constant term of ft by a subtraction of constant multiple of fi+\ for each 
0 < i < L - 2 and by a multiplication by w for i = L - \ : 


w 


1 -4 


1 


,yL-l 


1 

i fli-2 


/l-i 

fo 

fl 

fl-S 

/l-2 


= 0(w) = w 


fo 

f[ 


ifi 


L-U 


Eventually we obtain a new L-tuple /' = ^(/J,..., of power series from the original /. 

The above procedure is summarized as a left multiplication 


of an invertible matrix 


f = Tf 


0 

1 



1 


w 


a 

a 


1 

2 


1 


a 

a 


L-3 


1 


a 

a 


L-2 


(3.1) 


with det T = (-w)^“^. This is an analogue of the Euclidean algorithm and can be repeated gener- 
ically. Let f[k] = ^(fo[k],.. .,fL-i[k]) denote the corresponding vector of power series at the ^th 
step and let a'[k] denote their constant terms. Hence, we have 


f[k+l] = T[k]f[k] and /[0]=/, 


(3.2) 
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where 


0 

1 


Tm = 


1 

w 


a\k\ 

fl‘W 

1 


a\k\ 


On the other hand, solving (13.11) for / yields 


w 




(3.3) 


,/■=* 


(0 < / < L - 1) 


where = f^jw. Let us introduee the inhomogeneous coordinates (p = ^(^i,.. .,(Pl-\) by pi = 
fi! f{). Therefore, we have 


<Pi 

<Pl-i 


L-2 


-,L-\ 


+ W 


,/=' 


(0<i<L-2) 


(3.4) 


where (fo = 1. 


Definition 3.1 (cf. ifT^IIOll l. Let ip = ^(^i,..., Pl-\) be an (L - l)-vector such that p\ 4^ 0. Then 
the vector 

^ \^i (fil 

is called the reciprocal ofp. Note that = id. 

Under this notation, the correspondence (13.41) can be translated into 

t/ 1 4l-2' 


<Pl-\ <Pl-\ 


4l-\ 


= I ^{p) = a + wBp' 


or equivalently into 


where 


‘P 


1 


a + wBp' ’ 


^L-l 




n ^ ^ ... 1 

a' 4 

1 a* a* 

^ 4 ■ ■ ■ 

r 

a' 


a^-\ 

and B = 

1 / _7 

aL L 



1 


Namely, the vector a is determined as the constant term of l ^{p) and the matrix B is then specified 
by a. Let p\k\ denote the inhomogeneous coordinates of the vector f\k\ 6 C^[[w]] at the kih. step. 
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Taking the reciprocal repeatedly in this way, we obtain formally the vector continued fraction 




a[0] + 


w5[0] 


«[!] + 


w5[l] 


a[2] + 


w5[2] 
a[3] + --. 


(3.5) 


which is regarded as an (L - l)-dimensional generalization of the Stieltjes-type continued fraction. 
Refer to IfTTl Appendix A] for a classification of continued fractions. Our algorithm differs from 
the other known examples such as the Jacobi-Perron algorithm; cf. [|20l|2Tl. Note also that some 
dynamical system, like the Toda lattice, has been studied based on the connection among the 
Jacobi-Perron algorithm, rational approximations and bi-orthogonal polynomials; see [fOU . 

The following theorem can be verified straightforwardly through the above algorithm, as well 
as the case of a Stieltjes-type continued fraction (i.e., L = 2 case). 


Theorem 3.2. The kth convergents {rational functions) 


n, 


a[0]’ 


n. 


6 &-\w) 


a[0] + 


wfi[0] 


«[!] + 


wfi[l] 


<3([2] -l- ■. -l- 


wB\k - 2] 
a[k- 1] 


(k>2) 


of the vector continued fraction (13.51) provide approximants of the vector (p e Cf ^ [[w]] of power 
series in the sense that p -XVk = 0{w^). 

In calculating Hk, it is convenient to apply the projective transformations (13.31) successively as 
follows: 


TUq 


T 


UJ \ 

V7 1 

= r[0]“^r[i]"' - i]“' 

0 

and n^. = — 

VTq 

TZT2 



0 




3.2 Schlesinger transformations, revisited 

Let y = T(z) be a solution (11.81) of the Fuchsian system (11.71) having the local behaviors 

Y = 'T(z) • diag (z''°'Oo<i<L-i (near z = 0) 

= <I>(w) • diag (w'^°°’^)o<;<L-i ■ C (near z = 1/w = oo) 
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where the power series parts are normalized by 


'P(z) = 

1 • 

* 

- 1 - 0(z), 0(w) = 

* 




1 



* 


+ 0{w) 


and C is the conneetion matrix. Let / = ^(/o,..., /l-i) denote the first column of 0(w). 
It is clear from the construction of the matrix T = r[0] that 


ro 


/ 

* 






- 1 - 0{w) 

\ 


* 

/ 


diag(l,w \ 


near w = 0 (z = oo). On the other hand, it holds that 


T^ = z 


/ 

^ • • • * 1 

'V 


1 * • • • * 0 



1 : 

+ 0(Z) 


* 0 


V 

1 0 

/ 


diag(l,...,l,z ') 


VL 


1 • • • * 

1 


+ 0{z) 


0 

1 0 


1 0 


near z = 0. After we repeat the same procedure L times, the power series part thus recovers its 
original form: 

(\l ••• *1 

r[L-i]---r[i]r[0]'P = z^“^ 


VL 


1 


+ 0(z) 


In conclusion, the matrix z^~^T[L - 1] • ■ • r[l]r[0] turns out to be a polynomial in z and to be 
the multiplier of the Schlesinger transformation shifting the characteristic exponents at z = oo by 
(L - 1, -1,..., -1); cf. Theorem [L^ 

Remark 3.3. In fact, the same approximation problem considered in Sect. [H appears in the follow¬ 
ing manner. Concerning the polynomial matrix w^T[L- 1] • • • r[l]r[0], we observe from the form 
of T[k] that the diagonal entries are all monic linear functions, the strictly upper triangular part is 
linear and divisible by w, and the strictly lower triangular part is a constant. Moreover, in view of 
(13.21) we have 

w^T[L-l]---T[\]T[0]f = Oiw^), 
which coincides with the approximation condition (11.11) . and thus 


w^r[L - i]---r[i]r[0] 


~( 0 ) 

Q (w) 


Q (w) 


under n = 1. 
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4 Application to isomonodromic deformations 

In this section we first review some basic results on the Schlesinger system, which governs isomon¬ 
odromic deformations of a Fuchsian system of linear ordinary differential equations. As explained 
in Sect. I1.4[ a Schlesinger transformation preserves the monodromy of the Fuchsian system un¬ 
der consideration and, thereby, leads to a discrete symmetry of the associated Schlesinger system. 
Combining this fact with the result in Sect, [^reveals a determinantal nature of isomonodromic de¬ 
formations. Next we treat a particular case of the Schlesinger systems unifying various Painleve- 
type differential equations and show its relationship with certain hypergeometric functions, which 
will be needed later. 


4.1 Schlesinger systems and their symmetries 

Let us consider again the Lx L Fuchsian system (11.71) : 


— 

dz 



(4.1) 


where Mq = 1 and = 0. We start with a well-known result on isomonodromic deformations of 

dni). 


Theorem 4.1. The monodromy of a fundamental solution Y, i.e. det Y 0, does not depend on 
u = (mi, ..., M/v) if and only if 

dY 

Bi = —(1 < / < N) (4.2) 


are rational functions in z. 


We henceforth impose on our Fuchsian system (14.11) the following assumptions: 


(i) all the residue matrices A, are semi-simple, i.e. diagonalizable; 

(ii) there is no integer difference other than zero among the eigenvalues of each A, . 

Let us choose a normalization as before such that A^+i and Aa ?+2 = - upper and lower 

triangular matrices, respectively. Then we can take a fundamental solution Y = Y{z) of the form 


Y = 'P(z) • diag (near z = 0) 

= 0(w) • diag (w''”Oo< 7 <l-i • C (near z = 1/w = oo) 


with 


ri ••• *1 


'F(z) = 


+ 0{z), 


1 



> 



1 


<l>(w) = 

^ ... ^ 

0(w) = 

1 

* • • • 1 

- 1 - 0{w) 

/ 
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and E = diag • • • > where eaeh ^ 0 may depend on « = (ui,..., u^)- Moreover, we 
have 

Y = Giil + 0(z - Ui))(z - Uit' (4.3) 

near each of the other regular singularities z = w, {0 < i < N), where G, and A, are certain constant 
matrices satisfying G,A,Gr^ = A,. Therefore, the monodromy matrices of Y attached to loops 
around z = m, (0 < i < N) and z = 0 , oo read 

^2;rV^ diag (£oj)o<j<L-l d\ag (Soo,j)o<j<L-l Q 


Suppose now that every monodromy matrix of Y is constant with respect to u and, additionally, 
so is the connection matrix C. Then the rational functions Bi = 5,(z) can be explicitly written as 


Ai 1 

Bi = — -(AOlt, 

Ui - z Ui 


(4.4) 


where (A,)lt denotes the lower triangular part of A,; see Appendix for details. Note in particular 
that the diagonal part (A, )d of A, is expressible in terms of E as 


The compatibility condition 


(AOd = -M;—logs. 

OUi 


^-^ + [A,Bi]=0 
OUi OZ 


(4.5) 


of (14.11) and (14.21) is equivalent to a set of nonlinear differential equations for the matrices A, with 
respect to u, which is called the Schlesinger system Il22l . If (L, A) = (2,1), then the Schlesinger 
system reduces to the sixth Painleve equation Pyi. 

Next we shall investigate how the solution Y and the coefficient A = A(z) of the Fuchsian 
system (14.11) are connected with each other. Concerning the power series expansion 


Y = 0(w) • diag (w®“'0 • C, 0(w) = ^ 


k=0 


at the point of infinity (z = 1 /w = oo), the coefficients turn out to be polynomials in the 


off-diagonal entries of the lower triangular matrix Ayv +2 = - Y!i=o A, through Frobenius’ method. 
Conversely, substituting this solution Y in (14.IL we find that 


dT 

“dT 


-w- 


dY 

) - 

dw 

d 


-w — (Oq -I- Oiw -I- diag (w'^‘”'0 ■ C 

dw ’ 

- [(Oo -I- Oiw -I- 02W^ diag (gooj) 

-1-0iw -I- 202W^ -I- 303W^ diag (w®”'0 • C 


and 


v+i 

i=0 


Ai 


N+l 


i=0 


iV+1 


UiW 

N 


-I- 


ZO 

i=0 


+ UiW + Ui^W^ -I- 


')Ai 


Y 
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Recall Un+\ = 0 here. Equating the coefficients of in these power series yields 

An+2 = ffio • diag (£oo,y) • Oo“' for k = 0. 


thus Aiv +2 becomes a polynomial in the entries of the leading coefficient Oq of the solution, and 
also 

N 

^ UjAi = (-Oi • diag {s^j + 1 ) + Oq ' diag (eooj) • Oq"' for = 1 . 

i=0 

Moreover, if one needs similar expressions for all other residue matrices At besides A^+i, it is 
convenient to use the deformation equation (14.21) : one can verify in fact 


A,= 


<90o 

dui 


Oo ^ 


Oi - 


dui 


Oo ^ 


(1 < / < N). 


In summary, each residue matrix A, of A{z) is expressible as a polynomial in the entries of the 
coefficients of Y (and their derivatives with respect to u), and vice versa. 

A Schlesinger transformation keeps the monodromy of the Fuchsian system invariant but shifts 
its characteristic exponents by integers; recall Sect. 11.41 Consequently, it gives rise to a discrete 
symmetry of the Schlesinger system via the above correspondence between the solutions and coef¬ 
ficients of the Fuchsian system. On the other hand, any ingredient of Schlesinger transformations 
or of the associated rational approximations is described in terms of block-Toeplitz determinants; 
recall Sect. l2l This fact thus provides a natural explanation for the determinantal structure ap¬ 
pearing in solutions of isomonodromic deformations, e.g. Painleve equations. Refer to [[81 for a 
detailed investigation of the determinantal structure in Jimbo-Miwa-Ueno’s T-functions IfTOl for a 
general framework admitting irregular singularities. 


4.2 Polynomial Hamiltonian system *Hl,n of isomonodromy type 

We turn now to a particular case of the Schlesinger systems, which will be the main subject in the 
rest of this paper. 

Consider an L x L Fuchsian system of the form (14.1!) whose spectral type is given by the 
partitions of L: 

1 , L - 1 at z = M, (0 < z < A) and 
1 , 1, at z = 0 , oo, 

L 

which indicate how the characteristic exponents overlap at each of the A -l- 3 singularities. Fix the 
characteristic exponents as listed in the following table (Riemann scheme): 


Singularity 

Characteristic exponents 

Ui (0 < i < A) 

(-0z,O,...,O) 

hn+i = 0 

{eo,ei,... ,eL-i) 

Un+2 = OO 

(Ko-eo,Ki -ei,...,/CL-l -fiL-l) 


Assume the sum of all the characteristic exponents equals zero (Fuchs relation), i.e. 


L-l N 


/=0 ;=0 


(4.7) 
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Let A^+i and A^+i = - be upper and lower triangular matrices, respectively. Then such a 

Fuchsian system, denoted by Xlw, can be parametrized as follows: 


Ai = ^ b[‘\ bf_^ ■ (c[|\ cf,..., with c[j^ = 1 (0 < / < N), 

eo Woj • • • wo,L-i 


^V+l - 


ei 


WL-2,L-\ 

et-i 


(4.8) 


under the relations 

(tr A; =) - 0; = ^ ^ “ X 


L-l 


N 


N 


k=0 


i=0 


i=0 


the last two of which come from the triangularity of Aj^+ 2 . Also, we can and will normalize the 
characteristic exponents at z = 0 by 


trAA,+i 




L- 1 
2 


(4.9) 


without loss of generality. 

As shown in [|^ . the Schlesinger system governing isomonodromic deformations of Xl,v 
reduces to the multi-time Hamiltonian system 'NlX- 

/ l<ij<iv ] 

dxj dpf ’ dxj dqf \ i<k<L -I )' 

Here we let x, = 1 /w, and define the Hamiltonian function //, by 


XiHi = euqfpf + ^ qfp^^q\^^pf + 

k=0 j=Q 0<k<l<L-l j=Q 

j*> 



E (i) ij) (_/) (i) 

q^Pk ^i Pi 


with xo = qf^ = 4'^ = 1, pT = ki - Ym=\ qfpf ^nd p^^ = 0,- - Yak=\ qfpf- Therefore, Hi is a 
polynomial in the unknowns {canonical variables) 


c^^ 

«('■) = and n® 

^k ( 0 ) Pk 




'k ^k 


The number of the constant parameters 


/ 1 < / < A \ 

\ 1 <)t<L- 1 /■ 


(4.10) 


{e, K, 6) = (eo, • • •, eL_i, /Co,..., kl-uOq, ■■■,0n) 


(4.11) 


contained in is essentially 2L -l- A - 1 in view of (14.71) and (14.91) . For example, the case where 
L = 2 and any A > 1 coincides with the Gamier system in A-variables and, thereby, the first 
nontrivial case ^-( 2,1 does with the Hamiltonian form of Fyi- 

Remark 4.2. We have a priori known from their spectral type that the Fuchsian systems equipped 
with the Riemann scheme (14.61) constitute a 2A(L- l)-dimensional family. The coordinates of such 
a family are called accessory parameters, which are realized by the 2A(L - 1) canonical variables 
(14.101) in this instance. 
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4.3 Solution of in terms of hypergeometric function 

Although the phase spaee of ^l,n is a quite-eomplieated algebraie variety in general, there exists 
a family of solutions parametrized by a point in the projeetive spaee when the constants 

(e, K, 6) take certain special values. In fact, these solutions are written in terms of the hypergeomet¬ 
ric function 


(Q^l)|m| • • • {aL-l)\m\ifi\)mi ' ' ' (fiN)mN ~ ' ~ ^2) 

(ri)|m|---(rL-i)|M| mi!---m;v! ’ 

where \m\ = mi H-hm/^ and {a\ = r(a + k)/r(a). If (L,N) = (2,1), then (14.121) is exactly GauB’s 

hypergeometric function. 

To state the result precisely, we introduce the integral representation 



3^0 



U(t) dti • • -dfL-i 
W\=\{h-\ - ti) 


L-l 


n 


r(a/)r(y/ - ai) 
r(r/) 


X ^L,V 



of Fl^n, where the multi-valued function U = U(t) in t = (t\, t 2 , ■■■, tt-i) is given by 


L-l N 

Uit) = Y]tr^‘^fti_,-tir-‘^'Y]^l-XitL^^^^^^ and to = 71 = 1 , 

/=i 1=1 

and the cycle c is chosen to be an (L - l)-simplex 


{0 < tL-l < ■ • ■ < t 2 < tl < 1} c 


(4.13) 


(4.14) 


Also, we introduce supplementarily the integrals 

(0 _ r U(t) dti • • -dtL-i / 1 < / < At \ 

“ 1 (xitz.-i -1) -ti) \ 1 < -(< A -1 / ■ 

m 

We are now ready to state the hypergeometric solution of FIl,n\ see [1251 Theorem 3.2]. 
Theorem 4.3. If kq - = 0 then the Hamiltonian system Flip possesses a solution 

q'^ = 0 and pf = 0,— 

7o 


under the correspondence 


tTk — ^0, pi — hi, yii — eii Co l^k 


of constant parameters. 

The vector-valued function 


J = 



( 1 ) ( 2 ) 
>7l^,7i , 


( 2 ) 


(V) 

,/l , 
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satisfies a certain linear PfafRan system of rank N{L - 1) + 1, whose fundamental solution is 
prepared by collecting admissible cycles along with the foregoing (L- l)-simplex (14.141) . Note that 
the linear space of these cycles, i.e. twisted de Rham homology group, is generated by the chambers 
framed by the real section of branch locus of U = U(t); see [[I1|25]|. Of course, Theorem 14. 3 1 is 
valid for any solution of 'Pl,n- 

The Fuchsian system is specialized as kq - Oi = 0 and 




= Kkyo, 




bf = -0,. if = 09f, 

1 < / < \ 
\<k<L-l / 


(4.15) 


along the above hypergeometric solution of it thus becomes reducible. In fact, via the gauge 
transformation 


Y = 


N 


W 


Ao-eo 


]~[(1 - UiW)' 




i=l 


r. 


we have a solution of the form 


Y' = 


1 0 

/i 

fh-i 




/ 

1 






-I- 0{w) 

V 


■ 1 

y 


diag )o<;<i-i. 


where fk = fk{w) (f < k < L - f) are holomorphic functions at w = 0 dehned by the integrals 

r U{t) dti • • -dtL-i 

Jc (1 - wtt-i) - h) 

Hk 

Note that an (L - 1) x (L - 1) matrix W can be described by Thomae’s hypergeometric function 
For details we refer to Il27l . in which a curious coincidence between 'Pl,n+i and the Lax 
pair of 'Hi n, i.e. the pair of the original Fuchsian system (14.11) and its deformation equation (14.21) . 
is also discussed. 

Our aim here is to generalize Theorem 14. 31 by application of Schlesinger transformations start¬ 
ing from this hypergeometric solution at Ko-Yjf=i - 0- Notice that the Schlesinger transformation 
established in Theorem 1 1.61 shifts the constant parameters (14.1 II) as 


(/Co, /Cl,..., Kl-i) ^ (/Co + n{L - 1), /Ci - n,..., Kl-i - n) 

while all the others are unchanged; cf. the Riemann scheme (14.61) . Hence, by virtue of the algebraic 
relation between the solution of JLl,n and the canonical variables of "Hl^n (recall Sect. 14. II and also 
(14.81) and (14.101) ). we know in principle how to derive a solution {qf, pf) of at Kq - Xijii Oi = 
n{L - 1) for any positive integer n even though the resulting expression in this way will be terribly 
complicated. In the next section we explore this problem to achieve much simpler formulae for 
these special solutions. 
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5 Solutions of in terms of iterated hypergeometric inte¬ 
grals 

This section is concerned with the Schlesinger transform of the hypergeometric solution of 
We present its explicit formula by using the block-Toeplitz determinant whose entries are given 
by the hypergeometric functions. Key ingredients of the argument are the determinantal repre¬ 
sentations for the approximation polynomials; see Sect. [2l Moreover, we prove through Fubini’s 
theorem and the Vandermonde determinant that these block-Toeplitz determinants can be written 
in the form of iterated hypergeometric integrals. Our result will be summarized in Theorem l5.3l 
which is regarded as a generalization of Theorem 14.31 i.e. the previously known hypergeometric 
solution of ^l,n- 

5.1 Preliminaries 

Let /o(w), /i(w), ..., fi-iiw) be the functions defined by 


/o(w) = 1, 



If the cycle c is chosen such that < oo, then fk(w) is holomorphic at w = 0. For instance, it is 
enough to choose a bounded cycle as c. Accordingly, we have a power series expansion 


oo 


/^(w) = Yj = ho + h\w + h\w^ + • • • 


with the coefficients 



dtj • • • 


(5.1) 


c Y[\=i {U-i - ti) 


for 1 <k<L-\. Observe that each can be regarded as a moment 



pr(c) 


of the ‘measure’ 



upon the following notations: 


pr(c) = [ti-i I (L, ..., tL- 2 , A-i) 6 c), 
c\tL-\ = {(L, • • • , tL- 2 ) I {t\, ■ • • , ti-l, tt-i) 6 c}. 
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fk{w) 


Namely fk{w) (l<fc<L-l)is written as the Stieltjes transform 

dukis) 

/pr(c) 

of a funetion fXk = 

In parallel, we introduce the functions 


r d/uk 

JpT(c) ^ 


ws 


h 


also. We thus see that 


where 


d//o(tL-l) 

Lemma 5.1. The following linear relations {contiguity relations) hold: 


) _ r iL-dU{t) dti • • • dti_i 
’ -'c - h) 

= f s-’ d/uois), 

Jpr(c) 

'-'cltL-i 0/=! (h-1 h) I 


h)-Xih)^,=Uh% 


where ii denotes the down-shift operator with respect to fii defined by 

iiifii) = A' - 1 and biifij) = fij (i + j). 

Proof By definitions (15.11) and (15.21) it is immediate to verify these formulae. 


(5.2) 


(5.3) 

(5.4) 


□ 


Remark 5.2. If c is chosen to be the (L - l)-simplex (14.141) : 

{0 < ti-i < ■ ■ ■ < t2 < ti < 1} d 

both (15.11) and (15.21) are written by the hypergeometric function Let us introduce the function 


h = h 


a,/3 


:x 


L-I 

n 

/=! 


r(cr/)r(y/ - ai) 

r(7/) 


xF 


L,N 


a,p 

7 


Then it holds by definition that 


h) = h 


(Ti + 7 + 1,..., ak-\ + j + l,ak + j, ■ ■ ■, (Xl-i + i,P _ 


7i + 7 + 1,..., + 7 + 1,7^+1 + 7,..., 7 l-i + 7 

For example we have 






{Q<k<L-\). 


h) = h 


h) = h 


h] = h 




O'! +],■■■, (XL-\ + ffi 
Ti + L • • •, 7l-\ + 7 
ori + 7,..., ctl-i + ffi 
71+7+ 1,72 + 7, ...,7l-i + 7 

«! + 7 + 1, 0:2 + / • • • , «L-1 + ffi 
7i + 7 + 1,72 + 7 + 1,73 + 7, ■ • ■, 7 l-i + 7 






and so forth. 
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It is convenient to prepare the notation of the bloek-Toeplitz determinant (ef. Remark for 
any L-tuple of nonnegative integers: 


n = (no, Hi,..., ni_i) 6 (Z>o)^. 


Let |«| = 2^0^ n,. We set 


:= det 


(|«|,no) 


A‘l^_^{\n\,nk) hl - i ) 


= det 


k blocks 

Aj,|(«0, l«l) 


A\^^{nk-i,\n\) 


YAf-\{nL_M 


L — k blocks 


(5.5) 


for eaeh k (0 < k < L). If « = 0 = (0,..., 0), we fix A'^*^^(0) = 1. As well as (12.21) . the symbol 
Af(k, 1) denotes the kx I reetangular Toeplitz matrix for the sequenee whose top left eomer 

is h‘j, and hf = 0 if j < 0. The seeond equality in (15.51) can be verified easily from (12.31) . Heneeforth 
we suppose 

n = {0,n,... ,n) 

unless expressly stated otherwise. We will often abbreviate A®(«) for n = (0, n,..., n) as A®. 
Note that this convention is eonsistent with the deseription in Remark 12.2[ We also prepare the 
eanonical basis {eo,ei,.. .,eL-\} of Z^, i.e. 


Ck = ( 0 ,..., 0 , 1 , 0 ,..., 0 ). 

As seen in (14.151) . the coefficient 

v+i ^ 

Mz) = V -^ (Wo = 1, Un+Y = 0) 

to 

of the Fuehsian system Xl,v attached to the hypergeometrie solution of with kq - Gi = 0 
(see Theorem 14.31) is expressed as A,- = (0 < z < N), where 

=/jO .T(o,/ci,/c 2 ,...,/cl_i), = |l,-^,...,-^j, and 

= -e, • er\hl), ..., kr\K-')), = (i,o,...,o) for i < / < iv. (5.6) 

Cf. (14.81) . Our next task is applying the Sehlesinger transformation to this Fuehsian system. 
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5.2 Calculation of the Schlesinger transform (I) 

To derive the aetion of the Sehlesinger transformation, we need basieally to deal with (11.101) : 


v+i 


^ v~' 1 di? 1 

A ^ A = 2^ —^ = RAR-^ + —R~\ 


i=0 


Z - Ui 


dz 


Namely, sinee both R = R{z) and R ^ are polynomials in z, eaeh residue matrix A, ean be ealeulated 
by 

A, = R{ui)AiR-\ui) for 0 < i + 1. (5.7) 

However, thanks to (14.51) . it is rather easy to ealeulate the diagonal parts even in the general ease. 
First we will demonstrate it. 

The multiplier = R{z) of the Sehlesinger transformation (see Theorem 1 1.61) ean be written, a 
little more speeifieally than (11.91) . as 


R = w 


VL 


eS'^(O) 

* e?(0) 


* * 




+ 0(w) 


reeall (11.51) . Multiplying 


O = 


/ 

1 






+ 0{w) 

V 


1 

/ 


by R from the left yields R<t> = 6 • diag ”,..., w ”) with 


O = 


/ 

1 


\ 




+ 0(w) 

V 


1 

/ 


as shown in Theorem 1 1.61 where both S and S are diagonal matriees independent of w = 1/z. We 
mention, without fear of repetition, that the Hermite-Pade approximation eondition (11.11) assures 


Rf = 


PnL 

p\l 




^n(L 1 ) Qf higher order) 


with / = ^(/o,..., /l-i) denoting the first eolumn of O. We thus find the formula 


ss“i = diag (p°^, ^ 

/ A(i) Ad) A(2) A(i-l) 

= (-l) -diagH-l) 
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by virtue of RemarkCombining this with (14.51) under jc, = 1/m,: 





we arrive at the formulae 


d 

(A,)d - (A,)d = Xi— logdiag 

OXi 




Xi diag 


Ad)’A(^) 
£,.AW . A(i) D,A^i) • A(2) 


A^Ad) 


Ad)A(2) 


£),A(z-i) . A(^) 
A(r-i)A(2-) 


(5.8) 


for 1 < i < where D, denotes the Hirota differential with respeet to djdxt. 

Next we turn to the particular ease, i.e. the Fuehsian system Xlw upon the substitution (14.151) 
eorresponding to the hypergeometrie solution (see Theorem 14.31) . Write the residue matriees as 
A, = (0 < i < N). In order to reconstruet the eanonical variables ,pf) of the Schlesinger 

transform of it is only neeessary to know the quantities 

= to OS'SiV (5.9) 

beeause we have already known from (15.81) the diagonal entries (A,)^.^, = for 

1 < / < A; ef. (14.101) . Caleulations of (15.91) are a little eomplieated; therefore, we will separately 
earry out the eases i = 0 and z 0 in Seets. l53] and 15.4[ respeetively. The result ean be found in 
Seet.[531 


5.3 Calculation of the Schlesinger transform (II) 

Consider the row veetor c*^°^R“'(l) in view of (15.71) and Mq = 1, where = (1, -1 //Zq, ..., -1//Zq); 
reeall (15.61) . 


(i) The 0th eomponent It follows from the determinantal representation of the polynomial 
and Pf \w) = [/z^o°^]q for / > 0 (see Proposition [23] and its sequel) that 




Np(0) 


det 


^0 + ■ ■ ■ +^‘m-vK + ■ • ■ + ^L-2’ • • - >^0 + h[,h^^,0 
Al{n,m+ 1) 


A^ \n,m+ 1) 


Summation over I = l,2,...,L-lof this formula entails 


L-l 

y = J— 


det 


/jO _ /jO ^0 _ ^0 ^0 _ ^0 ^0 _ ^0 A 

AlXn,m+ 1 ) 
A^“'('2,m+1) 
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via the contiguity relation (15.31) . Hence the 0th component of ^(1) reads 

+ Cq) 

hilIsS^\n) 


1 h \ 1 

_V' =_det 


0 /=! 


■ ■ - ’"O 

Al(n,m+ 1) 


A^, ^{n,m + 1) 


Here we have used (12.71) . 


(5.10) 


(ii) The k{> 0)th component Similarly it holds for I > 0 that 


Pf\l) = 


NP 


(k) 


det 


hL + • • • + /z^ ,, hL + ■ • • + /z^ , 

0 m-1’ 0 m-2’ 

Al(n,m) 


A* \n,m) 
Al(n - I, m) 

li-i 




Taking a sum and using (15.31) thus yield 


L-l 


- 4 y = n ^ ro det 

Ktt 


. , /tQ ~ '^1 9 '^0 


/z^ /z^ 

Ai(zz,zn) 


Alin -l,m) 


Al;_\(n,m) 
(_l)nZ:+iA(^+i)(„ + go-ei) 
/z!!A®(n) 


(5.11) 


Here we have used (12.81) . 


The ratio of (15.101) and (15.111) leads to the formula 



A<^^(n)A(^+P(« + ep - eu) 
A(^)(«)A(°)(n + eo) 


5.4 Calculation of the Schlesinger transform (III) 

Consider for I < i < N the row vector c^^''R~^{Ui) in view of (15.71) and iz, = l/x,, which is nothing 
but the top row of the matrix i?“^(zz,) due to = (1,0,..., 0); recall (15.61) . 
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(i) The 0th component We are interested in the following determinant 


xrpfixi) = 


NP' 


i(0) 


det 


\,Xi,..., x”^ 

Al^(n,m+ 1 ) 


,A^-\n,m+l) 


Subtracting the jth column multiplied by x, from the (j + l)th column for j = m, m - 
sequentially and using the contiguity relation (15.41) . we thus obtain 


xrp^S\xd = 


NP' 


.( 0 ) 


det 


1 

Alin, 1) 

0,...,0 

^i(Al_An,m)\ 

xr'" 


Al_^{n,m)P 

Al;\n, 1) 

44...; 

= f- 

Np(0) ' 

det 

1 

Al;;\{n,m)l 


(_l)m^,-»Y,(A(0)(„)) 


(ii) The k{> 0)th component In the same way as above, we find 


xr^^pfixd = 


-m+l 


NP' 


ik) 


det 


\,Xi,...,Xi” 

Alin, m) 



Al_^in,m-l) 

Al;\n,m) 
Alin - l,m) 
A^:\in,m) 

y -m+1 

- ' f 

NpW ' 

det 

A’^-\in,m- 1) 
Al_^in- l,m- 1) 
All}^in,m- 1) 

Al;_\in,m) 


1 

_ Al-_\in,m-1) 





J/ 


A®(«) 

Hence we verify from the ratio of (15.121) and (15.131) that 

- Ck)) 


= (-ir.rr 


(A«»(«)) ■ 


5.5 Result 

Summarizing the above we are led to the following result. 


1,...,1 


(5.12) 


(5.13) 
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Theorem 5.3. Let n be a positive integer. If kq - ~ 1) ^hen the Hamiltonian system 

'Hl n possesses a solution (q^‘\p^‘^) given by 


-(0 


r^{i) ^(0 

^kPk 


-Xi 


-Xi- 


lsP\n + eo)€i - eu)) 

A(*;+i)(n + 00- euYi {M^\n)) ’ 

AW(n)A(^+i)(n) 


under the correspondence 


ttk — ^k Pi — bi, yii — ej( eo ki^ + n 

of constant parameters, where n = (0,n,... ,n) e (Z>o)^. 

Remark 5.4. We have in faet an alternative expression 


ftO) fS) — a Y- 
dk Pk ~ ^i^i' 




of the above solution with no use of the Hirota differentials. This formula ean be verified by 


ealculating in the same manner as Seets. l53l and l5.4t further details might be left to the reader. 


As seen in Theorem 15.31 we have eonstrueted a partieular solution of LIl,n expressed in terms 
of the block-Toeplitz determinant iSpfn) whose entries are given by the hypergeometric functions. 
Finally we shall rewrite A^^^(n) as an iterated hypergeometric integral: if we remember that h'. is 
defined as a moment (see Sect. 15. II) . then we observe that 

A®(n) = rdet ,..., , ■ • ■, 

J L l<i<no l<i<nk-\ l<i<nk l<j<nL-i 

L-1 n„ 

X n n 


fl=0 7=1 

^Jdet 


k 1 fla 

[■^OJ ]l<!<|«| j [■^IJ ]l<r<|n| 5 • • ■ ]l<;<|n| 1 ^ M ^a,j 

l<j<no i<j<ni l<j<nL-i\ 


L-1 n„ 


X n n ^biPaiSaj) 


(5.14) 


fl=0 7=1 

through Fubini’s theorem. The Vandermonde determinant shows that 


det 


m0,7 l<!<|n| ’ p^lj l<!<|n| j • • • 5 ‘^L-l.y l<i<|n| ~ I ^c,d)t 

l<7<no 1<7<«1 I<7<«L-lJ (a,b)Hc4) 


where (a, b) > (c, d) is defined to mean 

“a > c” or “a = c and b > J”. 
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Also it holds that 


^ Sgncr ]~^ ^ - Sa,c), 

cr£Sym(n„) j=\ b<c 

where Sym(na) denotes the symmetric group on {1,2,.. .,na}. Therefore, because the value of 
integral (15.141) is invariant under a permutation of the variables {Sa,b}i<b<n„ for each 0 < a < L - 1, 
we conclude that 

Jk-l 


=n ;77 / n n n 

a=0 \a=0 j=l ) \{a,b)>(c,d) 


^c,d) 


L-l 


a=0 


I I I ^a,c) I I j). 

^ i=l 


V b<c 


A Verification of (SM 

Usually, we often normali z e the Fuchsian system (14.11) so that ^^^+2 = - A, (the residue matrix 
at z = oo) is diagonal when considering its isomonodromic deformation; see e.g. IfTOl . But, in this 
paper, we adopt a different normalization treating the two points z = 0 and z = oo equally; i.e. we 
choose (the residue matrix at z = 0) and A ^+2 to be upper and lower triangular, respectively. 
Note that the latter normalization is more versatile in a general setting than the former. In this 
appendix, for a supplement to Sect. 14.11 we demonstrate how to determine the coefficient B, = B,(z) 
of the deformation equation (14.21) : 

dY 

of(l4Tl). 

Differentiating (14.31) with respect to uj tells us that 
dY 

• if i j then —7”' is holomorphic at z = m,; 

OUj 

dY Ai 

• if z = j then — T"' =-^—h (holomorphic at z = zz,). 

dui z - Ui 


Similarly, we observe that 


dUj du; 


0 * 
0 


+ 0{z) near z = 0 


(A.l) 


and 


dY. 




50 


dui 


dui 


+ D(w) near z = 1/w = oo. 


(A.2) 


Here we have used the assumption that the connection matrix C does not depend on zz,. Conse- 
dY , 

quently, Bi = —7 is a rational function matrix in z that has only a simple pole at z = zz, with 

dUi 

residue -A, and thus 

A,- 


B, 


+ Ki. 


Ui-Z 
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where Ki is a eonstant matrix. Substituting z = oo (w = 0) in (IA.2I) shows that Ki is a lower 
triangular matrix. Therefore, substituting z = 0 in (lA.ll) . we eonelude that Ki = -(A;)lt/m,. 
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